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Abstract 

We calculate 800 coefficients of the high-temperature expansion of the mag- 
netic susceptibility of Dyson's hierarchical model with a Landau-Ginzburg 
measure. Log-periodic corrections to the scaling laws appear as in the case of 
a Ising measure. The period of oscillation appers to be a universal quantity 
given in good approximation by the logarithm of the largest eigenvalue of the 
linearized RG transformation, in agreement with a possibility suggested by 
K. Wilson and developed by Niemeijer and van Leeuwen. We estimate 7 to 
be 1.300 (with a systematic error of the order of 0.002) in good agreement 
with the results obtained with other methods such as the e-expansion. We 
briefly discuss the relationship between the oscillations and the zeros of the 
partition function near the critical point in the complex temperature plane. 
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I. INTRODUCTION 



A possible way of testing our understanding of second order phase transitions consists 
in calculating the critical exponents as accurately as possible. Ideally, one would like to 
use several independent methods and obtain an agreement within small errors. The renor- 
malization group method IE has provided several approximate methods to calculate the 
critical exponents of lattice models in various dimensions. On the other hand, the same 
exponents can be estimated from the analysis of high-temperature series ||. Showing that 
these methods give precisely the same answers has been a challenging problem M. In gen- 
eral, one would expect that a well-established discrepancy could either reveal new aspects 
of the critical behavior of the model considered or point out the inadequacy of some of the 
methods used. 

In order to carry through this program, one needs to overcome technical difficulties 
which are specific to the methods used. An important problem with the high-temperature 
expansion || is that one needs much longer series than the ones available |7],^J (which do 
not go beyond order 25 in most of the cases) in order to make precise estimates. On the 
other hand, a problem specific to the renormalization group method is that the practical 
implementation of the method usually requires projections into a manageable subset of 
parameters characterizing the interactions. 

It is nevertheless possible to design a non-trivial lattice model ||, referred to hereafter 
as Dyson's hierarchical model (in order to avoid confusion with other models also called 
"hierarchical"), which can be seen as an approximate version of nearest neighbor models 
and for which these two technical difficulties can be overcome. For Dyson's hierarchical 
model, the renormalization group transformation reduces to a simple integral equation in- 
volving only the local measure. This simplicity allows one to control rigorously || the 
renormalization group transformation and to obtain accurate estimates of the eigenvalues 



of the linearized renormalization group transformation [10|. More recently, we have shown 



that the recursion formula can be put in a form njl2| suitable to the calculation of the 
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high-temperature expansion to very large order. Consequently, Dyson's hierarchical model 
is well suited to compare the e-expansion and the high-temperature expansion. Note that 
unlike the e-expansion, the high-temperature expansion depends on the choice of a local 
measure of integration for the spin variables (e.g. a Ising or Landau- Ginzburg measure). In 
order to make this choice explicit when necessary, we will, for instance, speak of Dyson's 
hierarchical Ising model if we are using a Ising measure. 

In a recent publication [|l^] , we reported results concerning the high-temperature expan- 



sion of Dyson's hierarchical Ising model. We found clear evidence for oscillations in the 
quantity used to estimate the critical exponent 7, called the extrapolated slope (see section 
III). When using a log scale for the order in the high-temperature expansion, these oscilla- 
tions become regularly spaced. We provided two possible interpretations. The first is that 
the eigenvalues of the linearized renormalization group are complex. The second is that the 
eigenvalues stay real but that the constants appearing in the conventional parametrization 
of the magnetic susceptibility should be replaced by functions of j3 c — (3 invariant under 
the rescaling of j3 c — (3 by Aj., the largest eigenvalue of the linearized renormalization group 
transformation. Hereafter, we refer to this explanation as "the second possibility". This 
second possibility has been mentioned twice by K. Wilson JTJ and developed systematically 
by Niemeijer and van Leeuwen [13[ . 



In Ref. 0], we gave several arguments against the first possibility. A more convincing 



argument is given in section VII: explicit calculations of the first fourteen eigenvalues of the 
linearized renormalization group transformation not relying on the e or high-temperature 
expansion show no evidence for complex eigenvalues of the linearized transformation. In 
addition, all the results presented below support the second possibility. 

In this article, we report the results of calculations of the high-temperature expansion of 
the magnetic susceptibility of Dyson's hierarchical model up to order 800 with a Landau- 
Ginzburg measure. These calculations provide good evidence that the oscillations appear 
with a universal frequency given by the second possibility |I],|T3J discussed above but with 
a measure- dependent phase and amplitude. Before going into the technical details related 



to the analysis of the series, we would like to state additional conclusions. First, we found 
no significant discrepancy between the high-temperature expansion and the e expansion. 
Second, with the existing methods, the high-temperature expansion appears as a rather 
inefficient way to estimate the critical exponents of Dyson's hierarchical model. Third, the 
high-temperature expansion reveals small oscillatory corrections to the scaling laws which 
cannot be detected from the study of the linearized renormalization group transformation. 

These conclusions were reached after a rather lengthy analysis. The second possibility 
introduces potentially an infinite number of Fourier coefficients and it is useful to first work 
with simplified examples in order to develop a strategy to fit the data with as few unknown 
parameters as possible. Solvable models where the second possibility is realized were pro- 
posed in Ref. |15[ . These models are sometimes called "Ising hierarchical lattice models" 
and should not be confused with Dyson's models. Further analysis of these models shows 
that the zeros of the partition function in the complex temperature plane are distributed on 
the (very decorative) Julia set |]16[ of a rational transformation. In particular, it is possible 
to relate the oscillations with poles of the Mellin transform located away from the real axis 
at the ferromagnetic critical point. In addition, the calculation of the amplitude of oscilla- 
tion for these models illustrates a feature which we believe is rather general, namely that 
the oscillations tend to "hide" themselves: large frequencies imply (exponentially) small 
amplitudes. 

This paper is organized as follows. In section II, we specify the models considered and 
the methods used for the calculations. In section III, we explain how to estimate the critical 
exponent 7 using the so-called extrapolated slope 0. We discuss the effects of the new 
oscillatory terms on this quantity, using assumptions which are motivated in subsequent 
sections. In section IV, we show that despite a large amplification, the systematic and 
numerical errors on the coefficients play no role in our discussion of the oscillations of the 
extrapolated slope. This section also provides a test of our calculation method in an explicitly 
solvable case, namely Dyson's hierarchical Gaussian model. 

Inspired by the Ising hierarchical lattice models and the analytical form of the one- 
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loop Feynman diagrams for Dyson's hierarchical model, we designed a simple mathematical 
function with a singularity corrected by log-periodic oscillations. This function is defined 
in section V. Its power singularity, as well as the frequency, amplitudes and phases of oscil- 
lations, can be explicitly calculated. We then show that these quantities can be extracted 
in good approximation from a numerical analysis of the extrapolated slope associated with 
the Taylor expansion of the function about a non-singular point. In section VI, we apply 
the methods developed in section V to fit the extrapolated slope associated with the various 
high-temperature expansions calculated. The analysis is complicated by the fact that the 
1/m corrections to the large m expansion - m being the order in the high-temperature ex- 
pansion - are enhanced by a factor which is approximately 160. We start with 5 parameter 
fits, which give robust values for the critical exponent 7 and the frequency of oscillation u. 
From the study of the errors, one can design fits with one or two more parameters which 
have smaller systematic errors and which are reasonably stable under small changes in the 
fitting interval or in the initial guesses for the values of the parameters. 

The results of the numerically stable fits are discussed in section VII. The linear relation 
between uj and 7 predicted by the second possibility is well obeyed and the value of 7 is 
in good agreement with the value obtained with the e— expansion, which we have checked 
using independent methods. All results agree within errors of the order 0.002. We have 
thus succeeded in finding a theoretical framework in which the new and existing results fit 
together. Many questions remain: What is the origin of the oscillation? Can we calculate the 
amplitudes of oscillation directly? Are similar phenomena present for models with nearest 
neighbor interactions? If the example of the solvable Ising hierarchical lattice models can 
be used as a guide, these questions require a better understanding of the susceptibility in 
the complex temperature plane. These questions are briefly discussed in section VIII. In 
particular, we give preliminary results concerning the zeros of the partition function in the 
complex temperature plane which suggests an accumulation of zeros near the critical point. 
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II. RECURSIVE CALCULATION OF THE HT EXPANSION 



In this section, we describe Dyson's hierarchical model and the methods used to calculate 
the high-temperature expansion of the magnetic susceptibility. The models considered here 

have 2 n sites. Labeling the sites with n indices x n x±, each index being or 1 , we can 

write the hamiltonian as 

1 n r 

# = -o£(i) £ ( £ " • (2- 1 ) 

1 = 1 X„,...,X t + 1 Xl,....,Xl 

The free parameter c which controls the strength of the interactions is set equal to 2 1 ~ 2 / D 
in order to approximate a nearest neighbor model in D-dimensions. In this article we only 
consider the case D — 3. The spins <7( Xn ; . .. are integrated with a local measure which 
needs to be specified. In the following we consider the Ising measure, where the spins take 
only the values ±1, and measures where the spin variables are integrated with a weight 
e~ j4o " 2 ~' Bcr4 , which we call Landau-Ginzburg measures. In the particular case B = 0, we 
obtain a Gaussian measure. In the following we have used A = 1/2 with B — 0.1 or B — 1. 

The integrations can be performed iteratively using a recursion formula studied in Ref. 
. Our calculation uses the Fourier transform of this recursion formula with a rescaling 
of the spin variable appropriate to the study of the high-temperature fixed point [|llj. It 
amounts to the repeated use of the recursion formula 

Rt+i(k) = C wexp(-^) W ^)(^(^)) 2 , (2.2) 

which is expanded to the desired order in (3. 

The initial condition for the Ising measure chosen here is Rq = cos(k). For the Landau- 
Ginsburg measure, the coefficients in the k— expansion have been evaluated numerically. The 
constant Ci + i is adjusted in such a way that i^ + i(0) = 1. After repeating this procedure n 
times, we can extract the finite volume magnetic susceptibility Xn(/3) = 1 + b\ yn (3 + &2,n/5 2 + 
... from the Taylor expansion of R n (k), which reads 1 — (l/2)k 2 x n + •••• This method has 



been presented for the Ising measure in Ref. |TT1 and checked using results obtained with 
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conventional graphical methods [17]. In the calculations presented below, we have used 
ra=100, which corresponds to a number of sites larger than 10 30 . The errors associated with 
the finite volume are negligible compared to the errors associated with numerical round-offs 
as explained in section IV. 

III. THE EXTRAPOLATED SLOPE 

In order to estimate 7, we will use a quantity called the extrapolated slope || and denoted 
S m hereafter. The justification for this will be made clear after we recall its definition. 
First, we define r m = b m /b m -x, the ratio of two successive coefficients. We then define the 
normalized slope S m and the extrapolated slope S m as 

S m = -m(m - l)(r m -r m _i)/(mr m - (m - l)r TO _i) ; 

(3.1) 

S m = mS m - (m - l)S m -i . 



In the conventional description [III] of the renormalization group flow near a fixed point 
with only one eigenvalue Ai > 1, the magnetic susceptibility can be expressed as 

X = iPc- /?)" 7 (^o + A x <fi e - + ....) , (3.2) 

with A = |/n(A2)|//n(Ai) and A2 being the largest of the remaining eigenvalues. It is usually 
assumed that these eigenvalues are real. When this is the case, one finds || that 

S m = 7 - 1 + Bm- A + 0{m- 2 ) . (3.3) 

Remarkably, the 1/m corrections coming from analytic contributions have disappeared, jus- 
tifying the choice of this quantity. Instead of this monotonic behavior, oscillations with a 
logarithmically increasing period were observed in Ref. Eq.(3.3) was then used, allow- 
ing B and A to be complex and selecting the real part of the modified expression. This 
introduces two new parameters, and the parametrization of the extrapolated slope becomes: 



S m = 7 — 1 — aim a2 cos(ujln{m) + a 3 ) . (3.4) 

This parameterization allows one to obtain good quality fits provided that m is not too small. 
This parametrization is compatible with two interpretations. The first one is that the 



eigenvalues of the linearized renormalization group are complex. We have given |T2[ several 
general arguments against this possibility and an explicit calculation reported in section VI 
makes this possibility quite implausible. The second possibility p]J13]] we have considered 
is that the eigenvalues stay real but the constants Aq and A\ in Eq. (3.2) are replaced 
by functions of f3 c — (3 invariant under the rescaling of (3 C — (3 by a factor (Ai) , where I is 
any positive or negative integer. This invariance implies that these functions are periodic 
functions in log{(3 c — (3) with period log(Xi) and can be expanded in integer powers of 
(P c — (3) in ( A i> . Consequently, we have the Fourier expansion: 

MPc-P) = EM&-/?) I ^ T • (3.5) 

ZGZ 

At this point, we have no additional information about these Fourier coefficients and 
possible restrictive relations among them. In the solvable examples [|Tj| where the second 
possibility is realized, the Fourier coefficients decrease exponentially with the mode number 
[TEL namely \au\ oc e~ UUJ ^ where 



^ = lit) ( 3 - 6 ) 

and u is a positive constant expected to be of order 1 but usually difficult to calculate. If a 
similar suppression occurs in the problem considered here, a truncation of the sum over the 
Fourier mode should provide acceptable approximations (see section V for an example). 

If we consider the new parametrization of the susceptibility - with the constants replaced 
by sums over Fourier modes - we obtain a parametrization of the HT coefficients as a linear 
combination of terms of the form (/3 C — j3) z . The asymptotic (at large m) form of the 
coefficients is obtained from 

{f3 c -(3f = (3: f; (")(-l) m (^) m (3.7) 

m=0 \ m / P c 



and the asymptotic form 

-z-1 



z \ , _ Nm m 



z + z 2 2z + 9z 2 + 10z 3 + 3z 4 Qz 2 + 17 z 3 + 17 z 4 + 7 z 5 + ^ 

x(1 + ^r + + is^ + ->- (3i 



From this we obtain the following asymptotic form for the coefficients: 

b m = m 1 ' 1 E K im il "(l + (( 7 + iIuj) 2 - ( 7 + iZw))/2m + ....) 
lez 

+m 7 - A - x E + ((7 - A + ilujf - (7 - A + iluj))/2m + ....) + ..., (3.9) 

where the i^; and L t are (unknown) coefficients proportional to the (unknown) Fourier 
coefficients. In the following, we consider the case of truncated Fourier series where only 
K , K±i and L are non-zero. A tedious calculation shows that to first order in Ki/K , L 
and 1/m, and neglecting terms of order L /m, we obtain 

S m = 7 - 1 + 2Re[i(u + uj 3 )m iuJ K t / K ] 

+L m- A (A 3 - A)2i?e[(A - A 3 - iuj + 3iA 2 uj + 3Auj 2 - i^m^i/lfo] 
+m- 1 i?e[(cj 2 + hiuj 2 - 2i^uj 3 + 7c/ - 2^uj 4 - ic^m^Xi/^o] . (3.10) 

From the solvable examples, we expect that \K 2 /K \ should be of the same order as \Ki/K \ 2 . 
The corrections of this order to S m read 



2Re[(2icu + 8iuj 3 )K 2 /K + (Aiuj 3 - iuj)^/ K ) 2 ] . (3.11) 

These corrections can be important at moderate uj (see section V). Importantly, we see that 
the 1/m terms have reappeared. In the case where uj » 1, we see that all the oscillating 
terms are approximately in phase and proportional to Re[im lw K\j K }. In the large uj limit, 
the 1/m corrections are enhanced by a factor oj 2 compared to the leading oscillating term. 
This feature will play an important role in the discussion of section VI. 

Before discussing the fits of the numerical values of the extrapolated slope for the Ising 
and the Landau- Ginzburg cases, we will first show that the errors made in the numerical 
calculations do not play any significant role and then discuss the fitting strategy with a 
solvable example. 
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IV. THE EFFECT OF VOLUME AND ROUND-OFF ERRORS 



In this section, we discuss the errors made in the calculation of the coefficients and show 
that they have no relevant effect on the extrapolated slope for the discussion which follows. 
There are two sources of errors: the numerical round-offs and the finite number of sites. We 
claim that with 2 100 sites and D = 3, the finite volume effects are several order of magnitude 
smaller than the round-off errors. 

From Eq. (2.2), one sees that the leading volume dependence will decay like (c/2) ra . This 
observation can be substantiated by using exact results at finite volume |T7| for low order 
coefficients, or by displaying the values of higher order coefficients at successive iterations 
as in Figure 1 of Ref. fll|| . In both cases, we observe that the (c/2) n law works remarkably 
well. For the calculations presented here, we have used c = 2s (i.e. D = 3) and n = 100, 
which gives volume effects on the order of 10 -20 . 

On the other hand, the round-off errors are expected to grow like the square root of the 
number of arithmetical operations. In Ref. ||11|| , we estimated this number as approximately 
nm 2 for a calculation up to order m in the high-temperature expansion with 2 n sites. As- 
suming a typical round-off error in double precision of the order of 10 -17 and n = 100, we 
estimate that the error on the m-th coefficient will be of order m x 10 -16 (or more con- 
servatively, bounded by m x 10~ 15 ). We have verified this approximate law by calculating 
the coefficients using a rescaled temperature and undoing this rescaling after the calcula- 
tion. We chose the rescaling factor to be 0.8482. The rescaled critical temperature is then 
approximately 1. This prevents the appearance of small numbers in the calculation. If all 
the calculations could be performed exactly, we would obtain the same results as with the 
original method. However, for calculations with finite precision, the two calculations have 
independent round-off errors. The difference between the coefficients obtained with the two 
procedures is shown in Fig. 1 and is compatible with the approximate law. From this, we 
conclude that for m < 1000, the errors on the coefficients should not exceed 10~ 12 . 

We are now left with the task of estimating the effects that the errors on the b m have 
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on S m . In general, S m is a function of b m , b m -i, b m -2, and b m -z- The numerical values of 
the derivative of S m with respect with these four variables are shown in Fig. 2 for a Ising 
measure. There is clearly a large amplification factor. From our upper bound on the errors 
on the coefficients, we conclude that the errors on the S m should be less than 10" 4 and will 
not play any role in the following. 

We have found independent checks of our error estimates. First, the smoothness of the 
data for the S m rules out numerical fluctuations which would be visible on graphs. The size 
of the data for the calculations with a Landau- Ginzburg measure allows a visual resolution 
of the order between 1CT 3 and 10~ 4 . Second, we have calculated S m in the Gaussian case 
where non-zero results are of purely numerical origin. The results are displayed in Fig. 3 . It 
shows that the numerical fluctuations for the Gaussian hierarchical model are smaller than 
10~ 7 for m < 200. This small number indicates that our previous estimates are conservative. 

The calculation of the large m coefficients requires a lot of computing time. We found 
that using a truncation in the expansion in k at order 100 could cut the computer time 
by a factor of order 100 while having very small effects on the values of the coefficients. 
If we plot the differences between the values obtained with the truncated and the regular 
method we obtain a graph very similar to Fig. 3. For m < 400, the differences are less 
than 4 x 10~ 6 , which is compatible with the numerical errors discussed above. The data for 
Landau- Ginzburg presented here has been calculated with the truncated method. 



V. DEVELOPING FITTING METHODS WITH A SIMPLE EXAMPLE 

The form of the coefficients given in Eq. (3.5) involves an infinite number of parameters. 
In order to see how one can obtain reasonable approximations with a manageable number of 
unknown parameter, we will first consider a simple example. One of the simplest examples 
of a function with a singularity and a log-periodic behavior is given by 

oo nn 

GW = £iw (5 ' 1) 



n 



This example has been motivated by the calculations of Refs. |16] and the form of the analytic 
expressions corresponding to one-loop Feynman diagrams for the hierarchical model. For 
definiteness, we shall only consider the case where A and B are real and A > B > 1. 

Picking an arbitrary positive value xq and introducing a new variable (3 = 1 — —we 
obtain the "high-temperature expansion" : 



oo 



with coefficients 



G{x) = £ b m (3 m , (5.2) 

n=0 



V , A A ^— . (5.3) 

f-i fl + A n Xn ) m+1 



^ (1 + A"x ) m+1 

The critical value of (3 is 1 and is obtained by setting x = in its definition. 



Using the Mellin transform technique discussed in Refs. ||16|| , we can rewrite 

G{x) = G reg (x) + G sing {x) , (5.4) 

with 

oo 

Greg(x) = £(-1)^(1 - BA^ 1 (5.5) 
1=0 

and 

Gsing{%) = ^X a sinMa+ipu)) ' (^) 

p=— oo v v 



where we have used the notation 



« = 7TT (5-7) 



and 

The complex part of the exponents comes from the fact that the Mellin transform of G(x) 
has poles away from the real axis. Substituting {(3 C — (3)xq for x and considering a as a critical 
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exponent, the analogy with the original problem becomes clear. Neglecting the regular part 
in (5.4) and proceeding as in section 3, we obtain the asymptotic form of the coefficients as 
in (3.9), with 7 replaced by a, Li — and 

K i = w t^^i • (5-9) 

r(a + zul)szn(n(a + iujI)) 

For large |/|, the magnitude of the coefficients decreases like e~^ u} ^\l\^~ a . One sees that 
fast oscillations have small amplitudes and vice-versa. This makes the oscillations hard to 
observe. In order to get an idea of how to obtain suitable truncations of the expansion given 
in Eq. (3.5), we have selected the values A = 3, B = 10 and x = 1 and calculated the 
coefficients with the exact formula (5.3). The sums were truncated in such a way that the 
remainder would be less than 1CT 16 . We then started fitting the corresponding S m using 
Eq (3.1). We first used a truncation where the Fourier modes with \l\ > 2 and corrections 
of order were dropped. We treated a, uj and the complex number as unknown 
coefficients and determined their values by minimizing the sum of the square of the errors 
with Powell's method. This allowed us to determine the order of magnitude of uj and a. 
Plotting the difference between the best fit and the exact values versus the logarithm of m 
shows oscillations twice as rapid as the oscillations in the fit. In other words, we needed 
the I = ±2 terms. With these terms included and using the data for m > 30, we obtained 
uj = 2.727 and a = 0.4772, in agreement with the exact values given by Eq (5.7) and (5.8), 
with three significant digits. The data and the fit are shown in Fig. 4. In this simple 
example, we found that each correction taken into account improved the quality of the fits. 
This is related to the fact that uj takes a not too large value. As we now proceed to discuss, 
a substantially larger value of uj implies a rather more complicated situation. 

VI. FITTING THE EXTRAPOLATED SLOPE 

We now discuss the fits of the extrapolated slope for Dyson's hierarchical model. The data 
is shown in Fig. 5 for the various measures considered. From the equally spaced oscillations 



13 



in the ln(m) variable, one finds immediately that oj is approximately 18. According to the 
exponential suppression hypothesis, this large value makes plausible that only the Fourier 
modes with |/| < 1 should be kept. This simplification unfortunately has the counterpart 
that for large ui the 1/m expansion is effectively a u 2 /m expansion, as explained at the end 
of section III. 

To be more specific, the relative strength of the leading oscillations and their 1/m cor- 

2 

rections is approximately For uo = 18, this means that for m = 162 the leading term 
and the first corrections have the same weight. In the example considered in the previous 
section, the critical value was m = 4 and good quality fits in the asymptotic region required 
considering the data for values of m larger than about ten times this critical value - which 
represents dropping only 5 percent of the data. For the hierarchical model, our data is 
limited to 5 times the critical value. Consequently, we probably need about 2500 coefficients 
in order to get results as accurate as in the example of section V. Despite this remark, an 
unbiased parametrization of the form 

S m = 7 — 1 + aim~ a ' 2 cos(uln(m) + a 3 ) + a 4 m _a2 + +a 5 cos(culn(m) + a 6 ) 

+a^mT l cos(uiln(m) + a 8 ) (6.1) 

gives very good quality fits provided that we disregard the low m data (see below). An 
example of such a fit is displayed in Fig. 6. The difference between the data and the fit is 
barely visible for m > 100. For m < 100 - where we do not have any reason to believe in 
the validity of the 1/m expansion - the frequency is still well fitted but not the amplitude. 
The assumption that only the Fourier modes with \l\ < 1 should be retained can be checked 
explicitly from the fact that the differences between the fit and the data do not show more 
rapid oscillations (unlike in the previous section, where the \l\ = 2 modes were important). 

If we want to have any chance of using the 1/m expansion as a guide, it is clear that 
we have to retain the data for m > m m i n , with m min larger than say 200. Varying m m i n 
and the initial values of the parameters provides a stability test. It appears that for the 10 
parameter fits mentioned above or their restriction to the 8 parameter case where all the 
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phases of the oscillatory terms are taken equal, the values of the fitted parameters depend 
sensitively on the value of m min and on the initial values. In particular, it makes no sense 
to check if the independent parameters satisfy relations dictated by the analysis of section 
3. 

We have nevertheless been able to design a stable procedure with less parameters. To 
assess the stability, we vary m min between 200 and 400, keeping m max at 800. The upper 
value of m min is chosen in such a way that we have at least two complete oscillations. We 
first set 04, 05 and a 7 equal to zero, which yields a parametrization of S m as in Eq. (3.4). 
These restricted fits do not suffer from the sensitive dependence mentioned above. We then 
analyze the errors as a function of to. In all the cases considered, the difference between the 
fit and the data is much smaller than the amplitude of the oscillations (for to > 200) and 
can be approximated by a constant plus a negative power of to. Putting together the fit 
of the extrapolated slope and the fit of the differences, we were able to obtain 6 parameter 
fits with a good stability and small systematic errors in 7. We now discuss the two cases 
separately. 

In the Ising case, the decay of the oscillations controlled by mr a ' 2 in the 5 parameter fit 
and the decay of the errors are both approximately m~ - 6 . We thus decided to use Eq.(6.1) 
with a 5 and aj equal to zero (making a 6 and as irrelevant). The 6 parameter fits so obtained 
are then reasonably stable under small changes in m min (see Figs. 8 and 9). Nevertheless, a 
systematic tendency can be observed: when m min is varied between 200 and 400, a 2 evolves 
slowly from 0.67 to 0.57. It is conceivable that if we had data at larger to, a 2 would evolve 
toward its expected value 0.46. 

In the Landau-Ginzburg case, the value of a 2 obtained from the 5 parameter fits is very 
small and the amplitude is in first approximation constant. We thus set a\ and a 7 equal to 
zero while a 5 parametrizes the amplitude of the oscillations and a 4 corrects the systematic 
errors. The power a 2 does not have the smooth behavior under a change of m min it had in 
the Ising case, however it does the job that it is required to do: the errors are small and do 
not show any kind or tilt or period doubling. These errors are displayed on Fig. 7. Their 
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order of magnitude is 10~ 3 , which can be used as a rough estimate of our systematic errors. 
Statistical errors due to the round-off errors are visible on the right part of the graph and 
clearly smaller by at least one order of magnitude. We now proceed to discuss the estimation 
of the most important quantities (7 and uj) from these fits. 



VII. ESTIMATION OF 7 AND w AND COMPARISON WITH EXISTING 

RESULTS 

The values of 7 as a function of m m i n are displayed in Fig. 8. The mean values are 1.3023 
in the Ising case and 1.2998 (1.2978) in the Landau-Ginzburg case with B — 1 (B — 0.1). 
We conclude that 7 = 1.300 with a systematic error of the order of 0.002. As explained in 
the previous section, a precise estimation of the subleading exponents seems difficult. 

Our results can be compared with those obtained from the e-expansion [1(J, namely 
Ai = 1.427 and A 2 = 0.85. These results imply 7 = 1.300 and A = 0.46. We have 
checked these results with methods which do not rely on the e-expansion or expansions in 
the renormalized coupling constants. First, we have adapted a numerical method discussed 
in Refs. [[[],[§ to the case of the hierarchical model. We obtained Ai = 1.426. Second, we 



have used a truncated and rescaled |TT| version of Eq. (2.2) which corresponds to the usual 
renormalization group transformation. Using fixed values of beta and retaining only terms 
of order up to k 28 at each step of the calculation, we were able to determine the fixed point 
and the linearized renormalization group transformation in this 14 dimensional subspace. 
Diagonalizing this matrix, we found Ai = 1.426 and A2 = 0.853. The corresponding value of 
7 is 1.302. For both methods, the errors can be estimated by comparing the linearizations 
obtained for successive iterations near the fixed point. The order of magnitude of these 
errors is 0.001 in both cases. As a by-product, we also found that all the other eigenvalues 
were (robustly) real, ruling out the possibility of complex eigenvalues. 

We now consider the values of uj. A distinct signature of the "second possibility" dis- 
cussed in Refs. [0,|T^,|r^,|T^] is the relation 
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This relation is well-obeyed by the a-priori independent quantities used in the fits as shown 
in Fig. 9. 



VIII. OPEN QUESTIONS AND CONCLUSIONS 

We have thus succeeded in finding a theoretical framework in which the new and existing 
results appear compatible within errors of the order of 0.002. In addition, we also have a 
qualitative understanding of the behavior of the extrapolated slope in the low m region. 
Many questions remain to be answered. First, we would like to understand the origin of 
the oscillations. If the example of the solvable Ising hierarchical lattice models can be used 
as a guide, the oscillations are due to poles of the Mellin transform located away from 
the real axis. These poles are related to an accumulation of singularities at the critical 
point. We have tried to get an indication that a similar mechanism would be present for the 
models considered here. As a first step, we have calculated the expansion of the partition 



function about {3 = 1.179, a good estimate |T3 of the critical temperature. We have carried 



the expansion up to order 10 for 2 n sites with n = 6 to 12. The zeros are displayed in 
Fig. 10. It appears that the approximate half-circle on which they lay shrinks around the 
critical point when the volume increases. It is not clear that the polynomial expansion is 
a good approximation. This could in principle be checked by searching for the exact zeros. 
However, this is a much harder calculation because due to the existence of couplings of 
different strengths, the partition function cannot be written as a polynomial in a single 
variable of the form e vf3 . 

The existence of log-periodic corrections to a singular behavior seems to be a feature of 
hierarchically organized systems. Empirical observations of such a phenomena have been 
suggested as a possible way to predict the occurrence of earthquakes [|19|] and stock market 
crashes [[2tJ . Are similar phenomena present for translationally invariant models with nearest 
neighbor interactions? Using the longest series available |§ for a nearest neighbor model, 
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namely the two-dimensional Ising model on a square lattice we found no clear evidence for 
regular log-periodic oscillations comparable to those seen in Fig. 5. However, the situation 
is complicated by the existence an antiferromagnetic point at (3 = —/3 C . We used an Euler 
transformation as discussed in Ref. || to eliminate this problem and found no indications 
of oscillations having a period that increases with m. On the other hand, the zeros of the 



partition function in the complex temperature plane have been studied [[nfl extensively. The 
zeros appear on two circles in the tanh{j3) plane, one of them going through the ferromagnetic 
critical point. Thus it seems incorrect to conclude that any accumulation of singularities will 
create oscillations. Approximate calculations of the Mellin transform of the susceptibility of 
the two-dimensional Ising model could shed some light on this question. 

We also would like to be able to calculate the amplitudes of oscillation with a method in- 
dependent of the high-temperature expansion. As explained in the introduction, the study 
of the linearized renormalization group transformation does not provide any indications 
concerning the oscillations. Up to now the global properties of the flows are only accessible 
through numerical approaches. The results presented here should be seen as an encourage- 
ment to develop and test global approaches to the renormalization group flows. 

We thank V.G.J. Rodgers for his help regarding computational problems. One of us 
(Y.M.), stayed at the Aspen Center for Physics during the last stage of this work and 
benefited from illuminating conversations with Martin Block concerning the fits of section 
VI. 
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Figure Captions 

Fig. 1: Difference between the S m calculated with the two procedures explained in the text. 
The solid line is m x 10~ 16 

Fig. 2: Derivatives of S m with respect to b m , 6 m _i, & m -2 and b m ^ 3 in the Ising case. 
Fig. 3: S m in the Gaussian case. 

Fig. 4: S m for the example of section V, with A = 3, B = 10 and x = 1 and the fit 
described in the text. 

Fig. 5: S m for the Ising model (crosses) and the Landau-Ginzburg model with B — 1 
(circles) and B = 0.1. (squares) 

Fig. 6: S m for the Ising model and a 10 parameter fit. 

Fig. 7: Difference between S m for Landau-Ginzburg with B = 1 and the fit given by Eq.(6.1) 
with 7 = 1.30137, = 17.716, a x = a 7 = 0, a 5 = -0.01084, a 6 = 0.3367, a 4 = 0.917 and 
a 2 =1.0589. 

Fig. 8: 7 as a function of m min , with m min between 200 and 400 by steps of 5 for the Ising 
model (circles) and the Landau-Ginzburg model with B = 1 (stars) and B = 0.1. (squares). 
Fig. 9: J"^) as a function of m min (as in Fig. 8) for the Ising model (circles) and the 
Landau- Ginzbug model with B = 1 (stars) and B = 0.1. (squares). 

Fig. 10: The zeros of the partition function in the complex temperature plane, in the Ising 
case with from 2 6 to 2 12 sites. The origin on the graph represents the point (3 = 1.179. The 
outer set of point (on a an approximate ellipse) is n = 6, the next set n = 7 etc.. 
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